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Abstract. We have obtained a high-S/N (900-1100), high-resolving-power (R=95000) spectrum of the metal-poor sub- 
giant HD 140283 in an effort to measure its 6 Li/ 7 Li isotope ratio. From a 1-D atmospheric analysis, we find a value con- 
sistent with zero, s Li/ 7 Li = 0.001, with an upper limit of 6 Li/ 7 Li<0.026. This measurement supersedes an earlier detection 
(0.040±0.015(lcr)) by one of the authors. HD 140283 provides no support for the suggestion that Population II stars may pre- 
serve their 6 Li on the portion of the subgiant branch where 7 Li is preserved. However, this star does not defeat the suggestion 
either; being at the cool end of subgiant branch of the Spite plateau, it may be sufficiently cool that 6 Li depletion has already 
set in, or the star may be sufficiently metal poor that little Galactic production of 6 Li had occurred. Continued investigation of 
other subgiants is necessary to test the idea. We also consider the implications of the HD140283 upper limit in conjunction with 
other measurements for models of 6 Li production by cosmic rays from supemovae and structure formation shocks. 

u 
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1 . Introduction (Deliyannis 1990), but not necessarily in slow-mixing models 

where the less-fragile elements 7 Li and 9 Be are depleted in con- 
Although the Big Bang is believed to be the major producer of cert with one another (De liyannis et al. 1998). Consequently, 
the 7 Li seen in Population II (Pop. II) stars, it is not believed reliable measurements of the Galactic evolution of 6 Li as well 
to be a significant source of the lighter isotope, 6 Li. A range as 7 Li could have wide implications. Since the fusion mech- 
of possible sites exist for 6 Li including not only spallative and anism (St eigman & Walker 1992) produces 6 Li without co- 
fusion sources normally associated with supernova-accelerated prod ucing beryllium or boron, observations of Be (Boesgaard 
Galactic cosmic rays (Walker et al. 1985; Steigman & Walker et al 1999) and B (Duncan et aL 1997) do not adequately con- 
1992), but also stellar flares (Deliyannis & Malaney 1995) and strain u productioni especially at the earlier epochs where fu- 
possibly shocks produced by large-scale-structure formation sion was more important than spallation. Models of spallative 
(Suzuki & Inoue 2002). Measurements of the 6 Li/ 7 Li ratio in production of Be and B are also affected by long-running uncer- 
metal-poor stars could therefore provide important constraints tainties in me Galactic chemical evolution of the most relevant 
on Li production following the Big Bang. The isotope ratio can heavy nucleus , oxygen. Consequently, the yield of 6 Li/ 9 Be de- 
also consttain possible destruction of Li, since 6 Li is destroyed pends not only on me cosmic ray energy spe ctrum (which may 
in stars at a lower temperature than 7 Li. 6 Li is more suscepti- include a substantia i i ow en ergy component (LEC); Vangioni- 
ble to destruction than 7 Li in some Li-depleting mechanisms Flam ^ Leh oucq, & Casse 1994), confinement (Prantzos, Casse, 
& Vangioni-Flam 1993; Fields, Olive, & Schramm 1994), the 
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* Based on data collected at the Subaru Telescope, which is oper- 1997) and composition (Fields et al. 1994; Ramaty et al. 1997; 
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Vangioni-Flam et al. 1999), but also on the fast particle and 
interstellar medium abundances. Depending on the particle 
source abundance distribution, the production ratio of 6 Li/ 9 Be 
ranges over 3-100 (Ramaty, Kozlovsky, & Lingenfelter 1996). 
Consequently, it is necessary to constrain 6 Li evolution directly 
from measurements of the 6 Li/ 7 Li ratio. 

Unfortunately, the measurement of 6 Li in stellar spectra is 
very difficult. The 6707 A transition, the only Li feature strong 
enough to permit an attempt, is a fine-structure doublet, and 
the isotopic displacement of the 6 Li lines from the 7 Li lines 
is comparable to both the doublet separation and to the in- 
trinsic line width. The latter is determined primarily by ther- 
mal Doppler broadening in the hot stellar atmosphere and by 
poorly characterised non-thermal (turbulent) motions. These 
factors, combined with the low fraction (<10%) of 6 Li, mean 
that high-resolution, high-signal-to-noise spectra are required. 
Even then, stellar models (Deliyannis 1990), supported by ob- 
servational evidence (Smith, Lambert & Nissen 1993), indicate 
that 6 Li is depleted below detection levels in all but the hottest 
main-sequence, Pop. II stars. Finding Pop. II stars hot enough 
to preserve 6 Li and bright enough to yield high S/N at high 
spectral resolution has been difficult. Progress was made with 
3- to 4-metre telescopes during the last decade: Smith et al. 
(1993, 1998) find 6 Li/Li = 0.06 + 0.03 in HD 84937 and 6 Li/Li 
= 0.05 ± 0.03 in BD+26°3578; Hobbs & Thorburn (1994, 
1997) measure 6 Li/ 7 Li = 0.08 ± 0.04 in HD 84937; and Cayrel 
et al. (1999) obtain 6 Li/ 7 Li = 0.052 + 0.019 in HD 84937. The 
reality of the 6 Li detection in HD84937 is supported by an anal- 
ysis of its K I 7698 A line (Smith et al. 2001). Also important 
are the numerous non-detections in other stars reported in those 
works. Isotope ratios have also been measured for a number of 
metal-poor disc stars (Nissen et al. 1999). The advent of 8-m 
telescopes should improve the situation. 

Deliyannis (1990, his Fig. 7) noted that 6 Li survival might 
also be high in subgiants containing Spite-plateau Li abun- 
dances, at least down to ~5800 K, since these stars had high 
effective temperatures on the main sequence. Deliyannis & 
Ryan (2000) set out to test this directly, and on the basis of 
a spectrum analysis by S.G.R. claimed a detection of 6 Li/ 7 Li at 
0.040+0.014(lcr) in the subgiant HD 140283. We set out to ver- 
ify this observation and analysis using the new High Dispersion 
Spectrograph (HDS) on the Subaru 8.2-m telescope (Noguchi 
et al. 2002). In the following sections of this paper, we present 
data with a S/N around 1000 per pixel, and report our finding 
of a very low upper limit, 6 Li/ 7 Li < 0.018, which supersedes 
S.G.R.'s earlier analysis. 

2. Observation 

During commissioning of HDS, we sought a very high-S/N ob- 
servation of HD 140283 to verify the earlier 6 Li detection at 
the same time as proving the capabilities of the new spectro- 
graph. The spectrum was obtained over two nights, 22 and 29 
July 2001, from thirteen exposures totalling 82 minutes. This 
gave a S/N = 900-1100 per 0.018 A pixel around 6707 A at 
R = A/AA = 95000, using a slit width of 0.4 arcsec (0.2 mm). 
The seeing was typically 0.55-0.6 arcsec, so the slit was well 
illuminated. This allows us to estimate the instrumental profile 



from the comparison spectrum, which we do below. The spec- 
trum was wavelength calibrated using a ThAr spectrum, and 
gave typical RMS errors of 0.0015 A. 

3. Spectral analysis 

3.1. General approach 

Our procedure was to compare the observed data with syn- 
thetic spectra computed using code originating with Cottrell 
& Norris (1978) and a ID model atmosphere grid by R. A. 
Bell (1983, priv. comm.). A model with T eS = 5750 K, logg 
= 3.4, [Fe/H] = -2.5, and microturbulence f = 1.4 km s _1 was 
used, along with the Li line list of Smith et al. (1998). Spectra 
were synthesised for various values of 6 Li/ 7 Li and convolved 
with macroturbulent and instrumental profiles. In order to do 
this we first determined the instrumental profile and then found 
(model-dependent) constraints on macroturbulence using pro- 
cedures described in Sects. 3.2 and 3.3. 

To determine the best fit when comparing synthetic and ob- 
served spectra, ax 2 test was employed, where 



and Oi is the observed continuum-normalised flux, 5, is the 
synthesised flux, and <x, is the standard deviation of the ob- 
served points defining the continuum. From this we calculated 
the reduced^ 2 , x% defined as 



where v is the number of degrees of freedom (e.g. Smith et al. 
1998). The best fit from a set of synthesised spectra is the one 
that minimises^ 2 . 

Once the observed spectra were normalised using neigh- 
bouring continuum windows, four variables affect the compar- 
ison between spectra: the abundance A(X) of element X, the 
wavelength shift AA, the macroturbulence F, and the 6 Li/ 7 Li 
ratio. To determine these, several iterations are needed as they 
are dependent on each other at some level. 

The wavelength shift AA that is allowed for each line may 
have any of several possible origins: an error in the applied 
Doppler correction for the star's motion; errors in the wave- 
length calibration of the ThAr frame (RMS = 0.0015 A); an er- 
ror in the wavelength of the line listed in the spectrum-synthesis 
linelist, which will differ from line to line; and possible varia- 
tions in the motions of different elements, ionisation states, and 
excitation states in a dynamic, 3D, real stellar atmosphere. We 
return to these shifts below, in Sect. 3.3. 

3.2. Instrumental profile 

The instrumental profile was calculated from a ThAr hollow- 
cathode-lamp spectrum over the interval 6660-6730 A taken 
with the same instrumentation setup as the stellar exposures. 
Nineteen emission lines of various strengths were isolated and 
normalised to the same strength by scaling them to the height 
of a fitted Gaussian. (In doing this we are not claiming that the 
ThAr profiles are Gaussian, merely that a Gaussian provides 



W. Aoki et al.: Lithium isotope ratio in HD 140283 



3 



a useful reference profile for normalising their heights.) On 
closer examination it was noted that six had weak lines in their 
wings, significant enough to impact weakly on our estimate of 
the instrumental profile, and so were rejected from the analysis. 
The remaining 1 3 unblended lines were then overlaid and aver- 
aged to give the instrumental profile shown in Figure 1(a). (The 
instrumental profiles for several setups have been given in 
Noguchi et al. (2002). Our result is similar.) It is comfort- 
ing that the lines give essentially identical profiles despite their 
different intensities. From the fall-off of the flux away from the 
line core, we believe that the lines effectively end around +0.3 
A, but we chose to truncate our calculations at 0.2 A. We esti- 
mate that < 0.5% of the flux lies beyond +0.2 A. 

Although we use the numerical profile in our stellar anal- 
ysis, it is instructive to investigate analytic approximations. 
Standard deviations were calculated over the intervals AA = 
±0.1 and +0.2 A, the former giving cr = 0.030 A and a tol- 
erable fit to the line core, and the latter giving cr = 0.037 A 
(Figure 1(b)). The larger standard deviation associated with the 
wider interval arises because 1.5% of the flux within AA = +0.2 
falls between ±0. 1 and ±0.2, outside the 3cr interval, whereas 
a true Gaussian would have less than 0.3% beyond 3cr. This in- 
flates the standard deviation if the Gaussian fit is evaluated out 
to ±0.2 A. This shows that the instrumental profile has broader 
wings than a Gaussian. The effect of (optionally) assuming a 
Gaussian instrumental profile is addressed in the stellar analy- 
sis below. 

Figure 1(b) also shows the bisector of the average instru- 
mental profile to illustrate its symmetry. Without magnification 
one would say it was perfectly symmetric. After magnification 
by a factor of ten, a small s-shape can be seen through the pro- 
file, but for our purposes the profile can be regarded as sym- 
metric. 

3.3. Macroturbulence 

We sought to constrain the macroturbulence from stellar spec- 
tral lines other than Li, and convolved the synthetic spectrum 
with a Gaussian broadening function of FWHM T\ km s _1 to 
model this. Synthetic spectra of five Ca I and Fe I lines (Table 
were calculated iteratively to find values of A(Ca) (or A(Fe)), 
AA and Fi that minimised^. As can be seen from column (6) 
of the Table, the values of Fi vary from line to line, but do not 
appear to correlate with atomic number, wavelength, excitation 
potential or equivalent width. It was noted that the Ca I 67 17 A 
line has a weak line in its long-wavelength wing that could not 
have been modelled simultaneously to the same precision, po- 
tentially degrading our results. Ca I 6717 also required a differ- 
ent AA value, -8 mA compared to -18 ± 2 mA for the others. 1 

1 Our analysis of the Ca I line assumed a wavelength of 67 17.688 A. 
Dr Bonifacio has kindly pointed out an unpublished measurement of 
the wavelength by Rosberg & Johansson, cited by Smith et al. (1998), 
giving the value 6717.677 A, which is 1 1 mA less than we have used. 
Adopting the Rosberg & Johansson value would change the AA value 
in Table 1 to -19 mA, which would be consistent with the other val- 
ues. However, it would not alter the broadening values for the line, 
which also seem slightly anomalous, possibly due to the presence of a 
weak blend as noted in the text. 
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Fig. 1. (a) points: Normalised flux versus distance from line 
centre for 13 ThAr lines. The peak intensity and central wave- 
length of each line was determined from a Gaussian fit. solid 
line: numerical instrumental profile obtained as the unweighted 
average of the points, (b) errorbars: numerical instrumental 
profile and standard error of each value, solid curve: Gaussian 
fit to instrumental profile over ±0.1 A. dotted curve: Gaussian 
fit to instrumental profile over ±0.2 A. vertical line: line bisec- 
tor of the instrumental profile, filled dots: line bisector enlarged 
xlO in wavelength. 



It was therefore decided to average the other four lines to give 
the preferred macroturbulence, (Fi ) = 4.65 km s . The consis- 
tency of the AA values for the four lines (excluding Ca I 6717) 
argues for a systematic error in the velocity correction applied 
during the data reduction as the source of the shift. 
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We also investigated the result of adopting a radial- 
tangential macroturbulent broadening profile, adopting the £r 
= </r and A R = A T prescription of Gray (1992, Chapter 18). 
The Ca I and Fe I lines examined above constrained £ RT 
to 4.21+0.13 km s _1 , whereupon a value 4.20 km s was 
adopted. 

To test the sensitivity of the results to the commonly- 
used assumption of Gaussian profiles, the best fit for a single 
Gaussian broadening function to model both the instrumental 
profile and macroturbulence was also calculated. The FWHM 
of this composite Gaussian is denoted r c . A similar spread 
in values for T c as for T\ was observed. These values gave 
(F c ) = 5.80 km s . Since F c represents a convolution of the 
instrumental profile with the supposed Gaussian macroturbu- 
lence F2, column (8) in Tabled estimates the macroturbulent 
portion assuming the instrumental profile is a Gaussian profile 
with a standard deviation of 0.03A (see § 3.2). 

3.4. Lithium isotope ratio 

With the instrumental profile determined numerically and the 
macroturbulence Fi constrained from four Ca and Fe lines, the 
synthesised and observed spectrum could be compared at Li 
6707 A. The^ 2 statistic was calculated over the 35 pixels from 
6707.50 A to 6708.20 A, where the Li profile is distinguish- 
able from the noise in the continuum. The continuum level was 
estimated out of this wavelength range. The effects of the un- 
certainties in the macroturbulence and the continuum level on 
the final results are estimated in subsection 3.5. 

The macroturbulence and the continuum level are two of 
the five free parameters determined by the profile fitting to the 
Li line in the analysis by Cayrel et al. (1999). However, the 
macroturbulence and the continuum level are, at least in prin- 
ciple, able to be determined independently of the Li line analy- 
sis. In contrast, the 6 Li and 7 Li abundances must be determined 
from the Li line analysis. Though the wavelength calibration 
and doppler correction can be made from the detailed analy- 
sis of spectral lines other than the Li line, there is still some 
uncertainty of the wavelengths of the line transitions. For this 
reason, the wavelength zero-point is also adjusted by the Li line 
analysis in the present work. 

Calculations of xl were performed covering a reasonable 
range of isotope ratios, total Li abundances, and wavelength 
shifts. These reduce the number of degrees of freedom by three, 
to 32. It was found during early work using multiple sampling 
of the independent variable x, that plots of xl vs x could be 
fit very well by a quadratic function, and hence the minimum 
value of x 1 ,- an d its corresponding x value could be found after 
a few calculations. This yielded 6 Li/ 7 Li = 0.001 as the best fit 
(see Tableland corresponding^ 2 -plots in Figure 2). 

The procedure was repeated for the pure Gaussian instru- 
mental profile, giving essentially the same result, 6 Li/ 7 Li = 
0.002. It can be seen from Figure 2 that adopting the pure 
Gaussian profile for the combined macroturbulent and instru- 
mental profiles is in close agreement with the preferred result 
based on the numerical instrumental profile. The calculation 




Q LJ I I I I I I I I I I I I I I I I I I I I 

0.01 0.02 0.03 0.04 
6 Li/ 7 Li 

Fig. 2. Xr values for five spectrum synthesis models as a func- 
tion of input lithium isotope ratio. Points indicate the actual cal- 
culations, while curves give the fitted quadratic functions used 
to identify the minimum, solid curve: preferred solution using 
numerical instrumental profile and Gaussian macroturbulence 
Ti = 4.65 km s constrained from four Ca and Fe lines, dotted 
curves: solutions in which F\ was varied to +lcr from the mean, 
to 4.50 and 4.80 km s _1 . dashed curves: Alternative solutions 
in which (i) pure Gaussian functions were used for instrumen- 
tal and macroturbulent broadening (T c =5.80 km s _1 ) and (ii) 
the HDS instrumental profile and radial-tangential macroturbu- 
lence (£rt = 4.20 km s ) were assumed. 

using the numerical instrumental profile also gives a marginally 
smaller^ 2 value. 

To test whether the Ca and Fe lines gave a reasonable value 
of Fj to use on the lithium lines, we performed a series of itera- 
tions on A(Li), A/16707 and 6 Li/ 7 Li in which Fi was also allowed 
to vary. The xl minimum was found for Fi = 4.53 km s _1 , 
shown in Table |2] This lies within lcr of the preferred value, 
4.65 km s _1 . Allowing T\ to vary used up another degree of 
freedom, but nevertheless put x 1 ,- marginally below the value 
for when Y\ took the preferred value constrained by the Ca and 
Fe lines. 

The lower broadening value weakens the wings of the pro- 
file and hence requires a higher 6 Li/ 7 Li ratio to fit the spe- 
cific observation, but 6 Li/ 7 Li remained below 0.01. Even al- 
lowing AA&jffj and A(Li) to vary, we found an almost linear 
relation between the 6 Li/ 7 Li ratio inferred and the adopted F\, 
6 Li/ 7 Li = -0.038F! + 0.1777 over the interval 4.50 < Fi < 
4.80; probably this relation holds over a much wider range of 
macroturbulent values. 

Figure 3 compares the observational data to the best fitting 
model and related models having 6 Li/ 7 Li = 0.02 and 0.04. 

We note that the equivalent width of the Li absorption line 
determined by the best-fit synthetic spectrum is 47.8 mA. This 
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4.10 
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6717.69 
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5.05 


6.00 


4.80 


4.60 
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6494.99 


2.4 


26.7 


-18 


4.60 


5.73 


4.78 


4.20 


Fel 


6678.00 


2.7 


12.4 


-15 


4.55 


5.63 


4.66 


4.15 
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± s.d. excl. 6717 




-18±2 


4.65±0.13 


5.79±0.17 


4.92+0.24 


4.21±0.13 


adopted 










4.65 


5.80 




4.20 



Notes: 

Ti = FWHM of Gaussian macroturbulence adopting numerical instrumental profile. 
T c = FWHM of composite Gaussian modelling macroturbulence and instrumental profile. 
T 2 = FWHM of Gaussian macroturbulence assuming r c incorporates a Gaussian 
instrumental profile with cr = 0.030 A. 



Table 2. Isotope ratios inferred from different spectral models 



Inst. Macroturb. 


Macrot. 


6 Li/ 7 Li 


xi 


A(Li) 


SA 


profile constraint 


km s" 1 








mA 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Preferred solution 












HDS Ca&Fe r 1= 4.65 
Other explored solutions 
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-0.001 


0.74 
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HDS Li 


r,=4.53 


0.005 


0.71 


2.181 


-28.0 


HDS Ca & Fe 


<Trt=4.20 


-0.007 


0.83 


2.188 


-29.2 


Gauss. Ca & Fe 


r c =5.80 


0.002 


0.81 


2.183 


-28.3 



value agrees well with that derived by Ford et al. (2002) for 
the same spectrum (48.1 mA). We also note that, even though 
the Li abundances derived by the present analysis are given in 
Table 2, the values are dependent on the choice of effective 
temperatures, while the isotope ratio is insensitive to the atmo- 
spheric parameters. 

3.5. Uncertainties 

We first estimate the statistical errors inherent in the spectrum 
with S/N^1000. We have tried to quantify this in two ways. 

The first method is to consider the interpretation of the 
X 2 statistic. The probability that a x 2 value as large as that 
measured should occur by chance is tabulated in many statis- 
tics books. Our plot of x 2 vs 6 Li/ 7 Li for the preferred analy- 
sis method indicates that very low xl values were obtained at 
the best fitting isotope ratio, corresponding to high probabili- 
ties that random statistical fluctuations could produce^ 2 values 
this large. The random fluctuations modelled in the^- 2 formal- 
ism are, we recall, the cr values associated with the noise in 
the continuum, in our case 0.0011. Away from the minimum, 
X 2 is larger and the probability that random fluctuations could 
produce such large x 2 values decreases. For 32 degrees of free- 
dom, the probability falls to 84. 1 % at x 2 = 1 -325, and to 97.8% 
atx 2 = 1-633. These correspond to 6 Li/ 7 Li = 0.019 and 0.023, 
respectively, for the adopted value of T - 4.65 km s _1 . 



However, due to the rebinning of the original spectra that 
occurs during data reduction, photon errors in adjacent pix- 
els are not fully independent. As the x 2 test assumes indepen- 
dent errors in adjacent pixels, the calculated chi-squared value 
is only an approximation to the value which ought to be min- 
imised, as pointed out by Cayrel et al. (1999) and Bonifacio & 
Caffau (2003). We continue to use the^ 2 procedure and in par- 
ticular select as our preferred parameter set that minimises the 
calculated^- 2 value, but we are unable to associate a confidence 
interval to any particular chi-squared values that we calculate. 

The second method involves conducting a Monte Carlo test 
in which a synthetic spectrum corresponding closely to the ob- 
served profile is subjected to Gaussian errors having the same 
distribution (cr = 0.0011) as the real data. A series of noisy, 
synthesised spectra based on a single input model are then sub- 
jected to the same ^-fitting procedures as the real analysis, to 
determine whether the input parameters can be recovered and 
with what accuracy. We performed 300 simulations to check 
the size of likely errors. Note that we are primarily interested 
in the spread of the results rather than the central value. For a 
synthetic spectrum having A(Li) = 2.223, a wavelength error 

= 0.001 A, and 6 Li/ 7 Li = 0.014, we inferred the follow- 
ing parameters: A (Li) = 2.224+0.001, AA = 0.0012+0.0011, 
and 6 Li/ 7 Li = 0.01 1+0.004. The inferred isotope ratios deviate 
from the input value with an RMS of 0.0042. This test sug- 
gests that, where the macroturbulence is constrained from other 
spectral features, the 6 Li/ 7 Li ratio can be recovered from data 
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i i i I i i i I i i i I i i i I i 

6707.4 6707.6 6707.8 6708 6708.2 
wavelength A 

Fig. 3. (a) Comparison of (curves) synthetic and (filled cir- 
cles) observed Li-doublet spectra for the preferred, best-fitting 
model in Figure 2 having 6 Li/ 7 Li = 0.00 and also models at 0.02 
and 0.04. Note that the higher 6 Li/ 7 Li ratios demand lower 7 Li 
abundances, resulting in a weakening of the line core. The inset 
shows three equal-area broadening profiles for (solid curve) the 
HDS instrumental profile, (dashed curve) T\ = 4.65 km s -1 , and 
(dotted curve) £rj = 4.20 km s _1 . The vertical dotted lines in- 
dicate the wavelength region over which the Li^ 2 analysis was 
conducted, (b) Difference in normalised flux between synthetic 
and observed spectra. Light, medium and heavy lines corre- 
spond to 6 Li/ 7 Li = 0.00, 0.02 and 0.04 respectively. The hori- 
zontal dotted lines correspond to + lcr in the observed flux. 



with S/N=900 at an RMS deviation = 0.0042. If the distribu- 
tion of the results of the above simulation is assumed to be the 
Gaussian, the 3cr limit is 0.013 in the 6 Li/ 7 Li ratio. 

The Monte Carlo estimate of the uncertainty in 6 Li/ 7 Li pro- 
vides a smaller error than that inferred from the x 2 statistic. 
Since the error estimate from the x 2 approach contains a diffi- 
culty in our case as mentioned above, we adopt the value es- 
timated by the Monte Carlo simulation (0.013 in the 6 Li/ 7 Li 
ratio) as the statistical error. 

It is easy to calculate the effect of uncertain T values in 
the inferred 6 Li/ 7 Li ratio because, as noted above, there is a 
linear dependence of the isotope ratio on the assumed value. A 
characteristic uncertainty <x r = 0.15 km s _1 translates to 0.006 
in the isotope ratio. 

We also estimated the errors due to the uncertainty of con- 
tinuum normalization by changing the continuum level in the 
analysis. We found almost a linear correlation between the 
assumed continuum level and resulting 6 Li/ 7 Li ratio. The as- 



sumption of by 0.1% higher continuum results in by -0.010 
lower 6 Li/ 7 Li ratio. 

Combining the limit ( 6 Li/ 7 Li <0.013)by the statistical error 
in quadrature with errors due to the uncertainties of macrotur- 
bulence and continuum level, we infer the upper limit of 6 Li/ 7 Li 
to be 0.018. 

4. Discussion 

The very low 6 Li/ 7 Li ratio found for HD 140283 is well be- 
low the value 0.040+0.015 inferred by S.G.R. from an earlier 
analysis of poorer data (Deliyannis & Ryan 2000). S.G.R's pre- 
vious analysis differed in several respects: the S/N was lower 
(435), and the sampling coarser (0.026 A pix -1 ), for the same 
resolving power. The reduction procedure also differed: in that 
work, an effort was made to constrain the fitting parameters 
independently rather than through a highly iterative x 1 analy- 
sis as employed here. It is probable that the present procedure 
has resulted in better optimisation of the fitting parameters than 
was achieved in the previous work. Application of the current 
technique to the older spectrum gave 6 Li/ 7 Li = 0.01 1 when F c = 
5.90 km s _1 was constrained by the Ca 1 6162 A line, and 0.034 
when F c = 5.38 km s -1 was constrained by the Li I 6707 A fit. 
These large variations in the inferred isotope ratio emphasise 
the need for good constraints on the macroturbulence if a reli- 
able isotope fraction is to be obtained. 

The analysis presented above is based solely on ID model 
atmospheres. Analyses using 3D model atmospheres, which 
have begun to appear in recent years (e.g. Asplund et al. 1999), 
offer the hope of replacing the empirical micro- and macro- 
turbulent parameters with physically-based velocity structures. 
The importance of this is clear from the strong correlation be- 
tween the input (adopted) macroturbulent broadening and the 
output isotope ratio. It can only be hoped that 3D calculations 
will soon become commonly accessible to all stellar spectro- 
scopists. 

Deliyannis' suggestion that sufficiently warm subgiants 
might preserve their main-sequence 6 Li complement finds no 
support from the current observation, but nor is the sugges- 
tion necessarily defeated. At = 5750 K, HD 140283 is very 
close to the effective temperature at which even 7 Li is seen to 
be diluted in subgiants (Pilachowski, Sneden & Booth 1993), 
so it is possible that 6 Li has begun to be depleted in this object. 
It is also amongst the most metal-poor of the stars in which 6 Li 
has been sought, so Galactic production may have been lower 
at this epoch. (Bear in mind, nevertheless, that 6 Li production 
via cosmic ray fusion is not as metallicity-dependent as pro- 
duction via spallation.) For these reasons it is still appropriate 
to examine warmer and/or higher-metallicity subgiants in an 
effort to detect 6 Li in them. 

However, it is also possible that our upper limit reflects a 
genuinely low 6 Li abundance at production. On the premise 
that depletion has not been significant, and that the HD 140283 
result is representative of stars in its metallicity range, we con- 
sider below the implications for different production models of 
Pop. II 6 Li. Fig. 4 shows our data together with previously pub- 
lished results of 6 Li detections in HD 84937, BD+26°3578 and 
G271-162 (Smith et al. 1998; Nissen et al. 2000 and references 
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therein). Also plotted are tentative detections for three addi- 
tional stars from Asplund et al. (2001), CD-30°18140, G13-9 
and HD 1606 17, achieved through preliminary analysis of re- 
cent VLT/UVES observations. Note that although HD 140283 
was additionally reported in Asplund et al. (2001) as a possi- 
ble detection, their latest analysis for this star is consistent with 
our limit here. Taken at face value, comparison of our upper 
limit with the other detections may suggest a relatively steep 
increase of 6 Li/H with metallicity near [Fe/H] ^ -2.4, followed 
by a slower rise. Alternatively, the data set may indicate typical 
6 Li abundances that are a factor of 2-3 lower than the highest 
measured values, which may be consistent with previous up- 
per limits for a few other stars (Smith et al. 1998, Hobbs et al. 
1999). 

The most widely discussed models so far for Pop. II 6 Li 
synthesis are based on spallation and/or fusion reactions in- 
duced by cosmic rays originating from supernovae (SNe). 
While they can successfully explain the Be and B observed in 
Pop. II stars, accounting for the observed 6 Li is rather prob- 
lematic; they require either an implausibly high cosmic ray in- 
jection efficiency (Ramaty et al. 2000, Suzuki & Yoshii 2001), 
or the presence of an additional low energy cosmic ray compo- 
nent lacking observational support (Vangioni-Flam et al. 1999). 
Even if the typical values of 6 Li abundances turn out to be 
lower than the highest measurements by a factor of 2-3, the 
data would still significantly exceed conservative predictions of 
such models assuming standard SN CR energetics and spectra, 
drawn as a dashed line (D) in Fig. 4 (see Suzuki & Inoue 2002). 
In some secondary SN CR models, a steeply rising log( 6 Li/H)- 
[Fe/H] relation can occur depending on the uncertain O-Fe re- 
lation (Fields & Olive 1999), but this is highly unlikely at the 
observed abundance levels in this metallicity range. 

A very different scenario for 6 Li production has recently 
been put forth by Suzuki & Inoue (2002): nuclear reactions in- 
duced by cosmic rays accelerated at structure formation (SF) 
shocks, i.e. gravitational virialization shocks driven by the in- 
fall and merging of gas in sub-Galactic haloes during hierar- 
chical build-up of structure in the early Galaxy. Such shocks 
are inevitable consequences of the currently standard theory of 
structure formation in the universe. Estimates for the specific 
energy dissipated at the main SF shock accompanying the final 
major merger give 6sf — 0.4keV per particle, compared to that 
for early SNe at €sn ~ 0.15keV per particle. Thus SF shocks 
can be potentially more energetic than SNe at early Galactic 
epochs, and the associated CRs can explain the 6 Li observa- 
tions more naturally. Since such shocks do not eject freshly 
synthesized CNO nor Fe, a - a fusion is the dominant pro- 
duction channel at low metallicities, which can generate large 
amounts of 6 Li with little Be or B and no direct correlation with 
Fe. A unique evolutionary behavior can arise, whereby 6 Li in- 
creases quickly at low metallicity (reflecting the main epoch of 
Galactic SF), followed by a plateau or a slow rise, in marked 
contrast to SN CR models for which log ( 6 Li/H) vs. [Fe/H] can 
never be much flatter than linear. Shown in Fig. 4 are three 
possible model curves for different parameter values of t SF , the 
main epoch of Galactic SF relative to halo chemical evolution, 
tsf, the main duration of SF, and jsf, the spectral index of 
injected particles (see Suzuki & Inoue 2002 for more details). 




-3 -2 -1 

[Fe/H] 



Fig. 4. Current observational data for 6 Li/H vs. [Fe/H] in 
Pop. II stars, compared with different models. The square, 
crosses and triangles correspond to the Subaru HDS upper 
limit for HD 140283, previous positive detections, and tentative 
VLT/UVES detections of Asplund et al. (2001), respectively. 
Model curves show the mean abundance predictions for SN 
CRs only (D, dashed), and SN plus SF CRs (A,B,C). Labels 
correspond to the following sets of parameters for t$F [Gyr], 
t sf [Gyr] and y SF : A (0.22, 0.1, 3), B (0.32, 0.1, 3) and C 
(0.22, 0.1,2). 

The current data set including the HD140283 upper limit 
may be consistent with curve C, corresponding to strong shocks 
with hard CR spectra (expected in case the preshock gas is ef- 
ficiently cooled by radiation), or a total SF CR energy a fac- 
tor of ~ 2.7 below the estimates described above (implying a 
lower total halo mass after the merger). If the steep rise of 6 Li 
is real, curve B may be a better representation, where the main 
SF shock at the final major merger occurs near [Fe/H] — 2, 
possibly being compatible with some other lines of evidence 
(e.g. Chiba & Beers 2000). For conclusive tests of the SF shock 
picture, more detailed and less parameterized modelling utiliz- 
ing e.g numerical simulations of galaxy formation is necessary. 
A further crucial prediction of the scenario is correlations be- 
tween the 6 Li abundance and the kinematic properties of Pop. 
II stars, which may also provide interesting insight into how the 
Galaxy formed. 

The low upper limit for HD 140283 is an important addition 
to the observational database for 6 Li in Pop. II stars. However, 
it is still insufficient for discriminating between different pro- 
duction models, especially in view of the uncertain extent of 
depletion. To elucidate the true origin of 6 Li as well as the ef- 
fects of stellar depletion in Pop. II stars, it is essential to obtain 
high quality data for a larger sample of stars with a wide range 
of metallicities and higher surface temperatures, an important 
goal for instruments such as the Subaru HDS. 

5. Conclusions 

We have obtained a high-S/N (900-1100), high-resolving- 
power (R=95000) spectrum of the metal-poor subgiant 
HD 140283 using HDS on the 8.2 m Subaru telescope, in an 
effort to measure its 6 Li/ 7 Li isotope ratio. We computed syn- 
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thetic spectra using a 1-D model-atmosphere code, and used a 
X 2 test to find the best fitting parameters. We find an isotope 
ratio consistent with zero, 6 Li/ 7 Li = 0.001, with an upper limit 
of 6 Li/ 7 Li<0.018 estimated by a Monte Carlo analysis. 

This measurement supersedes an earlier detection 
(0.040+0.015(lcr)) by one of the authors (Deliyannis & Ryan 
2000). HD 140283 provides no support for the suggestion 
that Pop. II stars may preserve their 6 Li on the portion of the 
subgiant branch where 7 Li is preserved. However, this star 
does not defeat the suggestion either; being at the cool end 
of subgiant branch of the Spite plateau, it may be sufficiently 
cool that 6 Li depletion has already set in, or the star may be 
sufficiently metal poor that little Galactic production of 6 Li 
had occurred. Continued investigation of other subgiants is 
necessary to test this idea. 

We also consider the HD 140283 upper limit along with 
other measurements in the context of 6 Li production models, 
particularly the structure formation shock scenario, in which 
unique evolutionary trends can be expected without any direct 
relation to the abundances of Be, B, CNO or Fe. If HD 140283 
has not been depleted in 6 Li, the upper limit can impose in- 
teresting constraints on the epoch and efficiency of dissipative 
dynamical processes that occurred during the formation of the 
Galaxy. Further observations together with improved theoret- 
ical modeling should allow us to test the structure formation 
scenario more critically and quantitatively, and to assess the 
value of 6 Li as a dynamical probe of the early Galaxy. 
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